Numerical Response of Owls to the Dampening of Small Mammal Population Cycles in Latvia

Simple Summary This article demonstrates the dampening of small mammal population dynamics and describes the numerical response of owls in Latvia. Numerical response was measured by diet, breeding performance, and trends in population change of six owl species. The responses varied among owl species, ranging from increased food niche breadth in more plastic species to reduced breeding performance and decreasing population size in more specialized species. The eagle owl seems to depend on voles in the previous autumn via the carry-over effect as measured by reduced breeding performance. Species more specialized in breeding in mature forests showed greater population declines, since mature forests are vital for owl breeding, as well as hold higher vole densities. Abstract Strong numerical and functional responses of owls to voles in cyclic environments are well known. However, there is insufficient knowledge from the boreonemoral region in particular, with depleted populations of small mammals. In this study, we describe the dynamics of the small mammal population in Latvia from 1991 to 2016 and link them to owl population characteristics. We used food niche breadth, number of fledglings, and population trends to lay out the numerical response of six owl species to dampened small mammal population cycles. We found temporarily increasing food niche breadth in tawny and Ural owls. There were no other responses in the tawny owl, whereas the breeding performance of three forest specialist species—pygmy, Tengmalm’s, and Ural owls—corresponded to the vole crash years in Fennoscandia. Moreover, the populations of forest specialist owls decreased, and the change in the Ural owl population can be attributed to the depletion of small mammal populations. We found evidence of a carry-over effect in the eagle owl arising from a strong correlation of declining breeding performance with the small mammal abundance indices in the previous autumn. We conclude that dampening of the small mammal population cycles is an important covariate of the likely effects of habitat destruction that needs to be investigated further, with stronger responses in more specialized (to prey or habitat) species.


Introduction
Small mammals play an important role in various ecological processes. This role ranges from influence on a natural succession [1] through influence on plant and microorganism community composition and chemistry [2] to demographic processes of small mammal predators [3][4][5][6] and even population processes and the behavior of directly unrelated species [7][8][9]. Fluctuations in the availability of small mammals are known to affect the demography of several species of birds of prey, and the analysis of the breeding performance of the latter can reveal large-scale spatiotemporal patterns of population dynamics of the former [10]. In the boreal region, small mammal populations typically Life 2023, 13, 572 3 of 28 three to four days [52]. Volunteers partially repeated this monitoring in 2015 and 2016. This scheme was conducted at four sites, but not every site was monitored every year (Table 1). REVIEW 3 of (a) (b) (c) (d)

Small Mammal Monitoring
Monitoring of the relative abundance of small mammals was conducted with sn traps from 1991 to 2016. It consisted of two schemes: the first with 2 transects per site, w 1 in a forest and 1 in a meadow (1); and the second with 11 transects per site, with 1 i meadow and 10 in different forest habitats (2).
The first scheme was officially run from 1991 to 2011. A total of 100 snap traps p transect (approx. 5 m between traps) were applied in autumn (August-September) three to four days [52]. Volunteers partially repeated this monitoring in 2015 and 20 This scheme was conducted at four sites, but not every site was monitored every y ( Table 1).
The second scheme was conducted in autumns (August-September) during the riod of 2012-2016 by volunteers. In this scheme, 20-25 snap trap transects (approx. 5 between traps) were applied in four areas (Table 1), although not all the areas were m itored every year. Forest transects were stratified into 10 categories as follows (minimu rotation ages for dominant tree species in Latvia are provided in Table A1): • YP-young (clearcuts and stands <7 years old) stands on poor soils; • YF-young (clearcuts and stands <7 years old) stands on fertile soils; • MPU-medium-aged (between 8 years and 80% of rotation age) stands on poor so without drainage; • MFU-medium-aged (between 8 years and 80% of rotation age) stands on fertile so   The second scheme was conducted in autumns (August-September) during the period of 2012-2016 by volunteers. In this scheme, 20-25 snap trap transects (approx. 5 m between traps) were applied in four areas (Table 1), although not all the areas were monitored every year. Forest transects were stratified into 10 categories as follows (minimum rotation ages for dominant tree species in Latvia are provided in Table A1):
The locations of all the small mammal monitoring sites are shown in Figure 1b.

Owl Diet during the Breeding Season
Owl diet analysis was based on prey remains and pellets found in nests or on the ground near the nest. Only the material from a single breeding occasion was used, based on annual nest-box and cavity inspections (GLAPAS, AEGFUN, STRALU, and STRURA) or based on the assumption that pellets cannot survive for many months in open nests or on the ground in the case of ASIOTU. The material was collected in autumn or winter from nest boxes and cavities and during chick ringing from ASIOTU nests. When collecting material, all the soft contents of nest boxes and cavities were removed. In the case of ASIOTU, all the useable material was collected. The distribution of the owl diet sampling sites in Latvia is shown in Figure 1c.
Analysis of the prey remains and identification of the minimal number of individuals were conducted as described by Vrezec et al. [35]. Insects were assumed to represent 1 g biomass, and amphibians and reptiles were assumed to weigh 16 g based on the average body mass of 100 measured individuals during chick ringing in 2016. In birds, the reference size group (i.e., woodpigeon, song thrush, chaffinch, and chiffchaff) weight from the general literature [53] was utilized. Region-specific weight of mammals from our trapping data or the literature [28,54] was utilized. We assumed young hare Lepus sp. to weigh 350 g.

Owl Population Change Monitoring
Monitoring of owl population change was conducted with traditional (territory mapping with playback broadcasting) methods [55][56][57][58][59] in permanent sample areas from 1991 to 2020, as well as with fully standardized point counts (with playback broadcasting) from national Breeding Birds of Prey Monitoring [60] from 2015 to 2022. The spatial distribution of monitoring sites is shown in Figure 1d.

Owl Breeding Performance
We used the number of fledglings per successful nest as a nest-level descriptor of the breeding performance strongly related to food availability. We used four data sources: nest-box inspections (1), information reported by ringers to the Latvian Ringing Centre (2), reports in the Breeding Birds of Prey Monitoring and the previous Monitoring for Owls (3), and citizen scientist reports on nature observation platform dabasdati.lv (4). The first two sources, as well as the third source (partially), covered information from the nests during ringing just before the young fledge. The citizen science and monitoring databases (partially) covered information on the number of young soon after fledging.
Most of the information before 2010 (apart from monitoring data) did not contain exact coordinates to be employed in spatial filtering and removal of duplicate records. Therefore, we combined nest-box inspection, monitoring, and ringing databases based on location attributes (indicated by nest name given by the ringer, which was most often also the person performing monitoring) to remove duplicates. We used citizen science reports only if there was no other information on the species in the particular year in the particular spatial reference. We used the national 1 km projected coordinate grid (epsg: 3059) if coordinates were known or the reported municipality otherwise as a spatial reference.

Data Analysis
We used R software [61] for data analysis, using the 'tidyverse' package [62] for data processing and visualizations and 'sf' [63] for spatial data. We treated results with p-values ≤ 0.05 as statistically significant but also reported insignificant results with full test descriptions.

Small Mammal Monitoring
We standardized the small mammal trapping data to the number of individuals per 100 trap days for further analysis and filtered only for autumn counts due to low representation of spring data. We used a graphical representation of the standardized counts per sampling area and habitat to compare variability between sites and habitats. We used generalized linear mixed-effects modelling (GLMM) to compare the differences in standardized densities between sampling areas, age classes, and soil fertility classes (including meadow habitats as a separate class in the latter two). We created two main effect models per comparison: • Random intercept per transect and the comparable variable in the fixed part; • Random intercept per transect and the comparable variable and year as a factor in the fixed part.
We utilized the Poisson family of distribution with a logarithmic link function and selected the best model based on the lowest value of sample-size-corrected Akaike information criterion value (AICc) [64]. We applied marginal means contrasting with Tukey's p-value correction for a post hoc analysis of the comparable variable between groups. For mixed-effects modelling, we implemented R packages 'lme4 [65] 'emmeans' [66] for contrasting.
We found no differences in either peaks or depressions between the sample areas in graphical analysis or mixed-effects models. Therefore, we used information from all the areas to obtain the countrywide small mammal population change index with TRIM analysis implemented in R package 'rtrim' [67]. The baseline model in this tool can generally be written as: where µ ij is an expected count, α i is a site parameter for site i, and β j is a time-point parameter for year j (for full explanation, see [68]). We created a model (model = 3) for the number of individuals pooled across species (small mammals) and separately for voles of genera Microtus and Clethrionomys. We evaluated serial correlation and overdispersion and selected the best model based on the lowest AIC value. We used the multiplicative slope of imputed values reported by TRIM to describe the overall population trend. We utilized graphical evaluation and Pearson's correlation analysis to compare the dynamics of both vole groups using yearly indices produced by TRIM analysis. We evaluated the presence of cyclicity with autocorrelation function analysis.

Owl Diet during the Breeding Season
We utilized only samples with at least 5 mammal prey individuals to avoid the influence of some very small samples. We used the small mammal data grouped according to genus as in [69] and other prey groups as described in Section 2.1.2. to calculate Levin's niche breadth (FNB): Life 2023, 13, 572 6 of 28 where p i is the fraction of a given prey in the total consumed biomass [70].
In Supplementary Materials, Table S1, we provide a description of the abundance and weight of the annual diet per owl species. For the description, we used the total number and cumulative biomass of the prey group, their arithmetical mean values, and proportions with Wilson's 95% confidence intervals.
We utilized linear regression (LM) analysis to evaluate the overall temporal change in food niche breadth. In species with repeated samples from the same nest sites, we also created linear mixed-effects models (LMMs) with nest ID as a random intercept. In the case of singularity in the random part, we employed the result of LM [71,72], as there were too few replicates per nest to contribute to the model fit. Model coefficients in the fixed part did not differ in these models. We fitted simple regression models with niche breadth as a response and year as a regressor variable. Both types of models were parameterized for Gaussian residual distribution with the identity link function.
Then, we evaluated the effect of each previously described small mammal population index on niche breadth with the same LM and LMM approach. We fitted models with the independent variable of the year of breeding, admitting a possibly reduced effect as the small mammal yearly indices represent late summer to early autumn rather than late spring to early summer, when owl breeding occurs.
Additionally, we used the biomass proportion of bank voles and voles of genus Microtus for regression analysis, as well as estimated population indices and bank vole proportion, with the index of Microtus voles. In this set of analyses, we applied GLMMs with a random intercept for nest ID and generalized linear models (GLMs) if the identifier did not contribute to model performance due to the number of replicates being too low [71,72]. We used the binomial family of distribution with logistic link function to compare the proportion of prey in the diet with its annual population index value in nature the same year.
Due to the small number of samples, we only described the AEGFUN diet without further analysis.

Owl Population Change Monitoring
Owl population monitoring started in 1990, but not all species had reliable information since the beginning of the period. Due to its preference for large forest massifs, STRURA was covered only since 1993. Due to a lack of knowledge on the monitoring of GLAPAS, its population change can be analyzed only since 2004. Due to the low population size, we did not have reliable data on the population change of BUBBUB.
We combined the data from both owl monitoring schemes if sites had all the planned census activities-for sample areas, sufficient coverage marked by an observer and for point count sites, four standardized visits to each point every year. To analyze population change, we drew on TRIM as described in Section 2.2.1. We exploited only the sites with information from at least two years and comparable effort according to the prerequisites [73].
We further calculated yearly indices and overall population change, as described by multiplicative slope with standard error [68,74], covering the whole available data period for the species. Then, we calculated two different slopes with relatively cyclic small mammal populations and since cycles had vanished. We used 2004 as the threshold for this division because it is:

•
The approximate time since when the small mammal populations did not recover to previous peaks; • The approximate midpoint of small mammal monitoring; • The approximate midpoint of STRURA monitoring; • The beginning of GLAPAS monitoring.
To calculate "before" and "after" trends, we selected the necessary parts of yearly indices and conducted linear regression on ln-transformed indices. We defined index 1991-2004 as "before" and index 2004-2016 as "after". To obtain significance tests, we transformed the time to start with year 1 in each group to use in regression with interaction between time and period. We defined "before" as a reference level.

Owl Breeding Performance
We calculated the annual mean values of spatially cleaned results because in most cases, we were not able to match different breeding performance reports to an exact breeding location (nest or territory). To provide more generalizable information, we calculated bootstrapped 95% confidence intervals from 1000 bootstrap resamples and implemented them in visualization. We used the annual mean values to establish a temporal trend of overall change in breeding performance. We employed slightly different approaches for further processing because the amount of information varied among owl species.
For the two species (STRALU and STRURA) with most data available on a nearly annual basis, we compared the temporal trends before and after the dampening of small mammal cycles in 2004. We used Gaussian linear regression with the annual mean number of fledglings as a response and compared the trends between the periods as in population change analysis.
We used all the years available to correlate the annual mean breeding performance of STRALU and STRURA with small mammal population indices in the year of breeding and one year before to evaluate a possible carry-over effect. In the case of ASIOTU and BUBBUB, we used a reduced timeframe (from 2002 and 2001, respectively) to avoid possible artefacts due to irregular sampling and a small sample size. We used Spearman's rank correlation analysis due to some outliers and a slightly curved scatter plot.
We harnessed the R package 'Hmisc' [75] for bootstrapping and base R for correlation and regression analysis.

Small Mammal Monitoring
The number of small mammals per 100 trap days over time in different sample areas and habitats is illustrated in Figure 2a. As the figure depicts, the peaks and depressions matched well between areas over time, with only a slight stochasticity between habitats within the same areas. This was confirmed by GLMM analysis, showing no significant differences in marginal mean ratios of the number of small mammals per 100 trap days among sample areas when accounting for an individual transect in a particular year (Table A2). There were observable differences in the relative abundances of small mammals between habitats ( Figure 2b). GLMM analysis revealed that meadow habitats had significantly lower abundances than any forest age group but with no differences between age groups (Table A3). Comparison of fertility groups revealed that meadows had significantly lower abundance and that forests on fertile soils had significantly higher abundance after accounting for multiple comparisons (Table A4). In every comparison, GLMM including the hierarchical random intercept of transect in year and only the variable of interest in the fixed part was the best-performing model (with the lowest AICc values).
We pooled all the results to conduct a population change analysis because there were no important differences between the areas. All three models suggest statistically significantly declining populations (total number of small mammals: S = 0.9671 ± 0.0083, We pooled all the results to conduct a population change analysis because there were no important differences between the areas. All three models suggest statistically significantly declining populations (total number of small mammals: S = 0.9671 ± 0.0083, p = 0.0007; Microtus voles: S = 0.9306 ± 0.0167, p = 0.0005; bank voles (Clethrionomys glareolus): S = 0.9706 ± 0.0128, p = 0.0325). Yearly abundance indices are shown in Figure 3a-c.  We pooled all the results to conduct a population change analysis because there were no important differences between the areas. All three models suggest statistically significantly declining populations (total number of small mammals: S = 0.9671 ± 0.0083, p = 0.0007; Microtus voles: S = 0.9306 ± 0.0167, p = 0.0005; bank voles (Clethrionomys glareolus): S = 0.9706 ± 0.0128, p = 0.0325). Yearly abundance indices are shown in Figure 3a-c.

Owl Breeding Season Diet
In total, 164 STRALU samples from 86 different locations covering 23 years, 56 STRURA samples from 38 different locations covering 15 years, 24 ASIOTU samples from 21 different locations covering 9 years, 7 GLAPAS samples from 7 different locations covering 7 years, and 2 AEGFUN samples from 2 different locations covering 2 years ( Figure 4, Table S1) Life 2023, 13, 572 9 of 28 were analyzed. The description of annual food composition per species is presented in Supplementary Materials, Table S1. ure 3d).

Owl Breeding Season Diet
In total, 164 STRALU samples from 86 different locations covering 23 years, 56 STRURA samples from 38 different locations covering 15 years, 24 ASIOTU samples from 21 different locations covering 9 years, 7 GLAPAS samples from 7 different locations covering 7 years, and 2 AEGFUN samples from 2 different locations covering 2 years ( Figure  4, Table S1) were analyzed. The description of annual food composition per species is presented in Supplementary Materials, Table S1. The overall average food niche breadth (FNB) of STRALU was 5.125 (95% bootstrapped confidence interval (bCI), 4.867-5.423). FNB increased significantly (β = 0.0840 ± 0.0198; t (129.2401) = 4.249; p < 0.0001) from 1992 to 2016. There were notable differences between samples in any given year (Figure 5a), and inclusion of nest ID as a random effect provided some help in dispersion taming (LMM: R 2 conditional = 0.115, R 2 marginal = 0.101, ICC = 0.015), indicating some degree of territory-specific variability. Although the explained variances were low, we found a statistically significant negative effect of each of the small mammal population indices on FNB ( Table 2). The proportion of voles (both groups) in owl prey correlated positively with the vole abundance indices, but STRALU showed preference for Microtus voles, as their abundance index had a significant negative correlation with the bank vole proportion in prey (Table 3). On average, Microtus voles accounted for 15.55% of biomass, whereas bank voles and voles in total accounted for 5.51% and 31.65% of biomass, respectively. The overall average food niche breadth (FNB) of STRALU was 5.125 (95% bootstrapped confidence interval (bCI), 4.867-5.423). FNB increased significantly (β = 0.0840 ± 0.0198; t (129.2401) = 4.249; p < 0.0001) from 1992 to 2016. There were notable differences between samples in any given year (Figure 5a), and inclusion of nest ID as a random effect provided some help in dispersion taming (LMM: R 2 conditional = 0.115, R 2 marginal = 0.101, ICC = 0.015), indicating some degree of territory-specific variability. Although the explained variances were low, we found a statistically significant negative effect of each of the small mammal population indices on FNB ( Table 2). The proportion of voles (both groups) in owl prey correlated positively with the vole abundance indices, but STRALU showed preference for Microtus voles, as their abundance index had a significant negative correlation with the bank vole proportion in prey (Table 3). On average, Microtus voles accounted for 15.55% of biomass, whereas bank voles and voles in total accounted for 5.51% and 31.65% of biomass, respectively.   STRURA also showed large variability in the diet, as overall FNB was 4.485 (95% bCI, 4.201-4.758). We found a significant (β = 0.0499 ± 0.0227; t (54) = 2.194; p = 0.0325) increase in FNB from 1994 to 2016 ( Figure 5). The overall variability was lower than in STRALU, but no nest-specific intercepts were found to improve the model, and LM could explain only about 6.5% of the variance (R 2 adj. = 0.06486). We found no correlation between FNB and small mammal population indices ( Table 2). The species showed a strong preference for Microtus voles, the proportion of which in prey correlated positively with its abundance index, while higher abundance in nature led to a lower proportion of bank voles in  STRURA also showed large variability in the diet, as overall FNB was 4.485 (95% bCI, 4.201-4.758). We found a significant (β = 0.0499 ± 0.0227; t (54) = 2.194; p = 0.0325) increase in FNB from 1994 to 2016 ( Figure 5). The overall variability was lower than in STRALU, but no nest-specific intercepts were found to improve the model, and LM could explain only about 6.5% of the variance (R 2 adj. = 0.06486). We found no correlation between FNB and small mammal population indices ( Table 2). The species showed a strong preference for Microtus voles, the proportion of which in prey correlated positively with its abundance index, while higher abundance in nature led to a lower proportion of bank voles in prey (Table 3). We suggest preference as a reason for the negative correlation of the bank vole abundance index with its proportion in prey because both vole abundance indices were correlated (Table 3). On average, Microtus voles accounted for 15.10% of biomass, whereas bank voles and voles in total accounted for 6.07% and 31.78% of biomass, respectively.
The overall FNB of GLAPAS was 3.526 (95% bCI 2.355-4.756), showing a temporal increase ( Figure 5). However, this increase was not found to be statistically significant (LM: β = 0.0420 ± 0.1136, t (5) = 0.369, p = 0.727; R 2 adj. = −0.1681), probably due to high intersample variability and small overall sample size. We found no correlation between FNB and small mammal population indices ( Table 2). Results of prey proportion and relative abundance in nature were similar to those for STRURA, but the preference for Microtus voles was greater, serving as better explanator of the bank vole proportion in prey (Table 3). On average, Microtus voles accounted for 10.03% of biomass, whereas bank voles and voles in total accounted for 11.31% and 29.46% of biomass, respectively.
The two analyzed samples of AEGFUN had FNBs of 1.588 and 4.318, respectively (Table S1). On average, Microtus voles accounted for 62.89% of biomass, whereas bank voles and voles in total accounted for 13.81% and 76.69% of biomass, respectively.

Owl Population Change
Our results differed between species when comparing owl population changes before and after small mammal cycle depletion ( Figure 6, Table 4). STRALU, with an overall (1990-2021) stable population (S = 1.002 ± 0.005), exhibited no significant difference in population trends before and after depletion (Table 4). Although the population experienced a considerable depression during the period of 2010-2012, it has since recovered (Figure 6a).
We obtained similar results for ASIOTU, with an overall (1990-2021) stable population (S = 0.992 ± 0.010) and no significant difference between the periods (Table 4). However, a visual extension of the trend since mammal population depletion suggested a decline that might be obscured by fluctuating population (Figure 6c).
The results of STRURA were different; although overall (1993-2021), the population was classified as stable (S = 1.014 ± 0.012), there was a significant difference in trends (Table 4). This species had a strongly increasing population before 2004 and a steep decline since small mammal depletion (Figure 6b).
Population change information for GLAPAS was available only since the depletion, and its overall (2004-2021) negative population trend (S = 0.965 ± 0.017) was similar to that observed in 2004-2016, reflecting a significant decline (Figure 6e, Table 4).
The results for AEGFUN are interesting, as the overall (1990-2021) population had a steep decline (S = 0.934 ± 0.020) that fit with estimated yearly indices (Figure 6d). Nevertheless, the difference between slopes "before" and "after" depletion was significant (Table 4) and suggests a steeper decline during the pronounced small mammal dynamics than since the depletion of cycles. However, visually extending the trajectory of "after" revealed a pattern to that in the "before" period; thus, the difference could be an artefact due to an increased influence of some years.

Owl Population Change
Our results differed between species when comparing owl population changes before and after small mammal cycle depletion ( Figure 6, Table 4). STRALU, with an overall (1990-2021) stable population (S = 1.002 ± 0.005), exhibited no significant difference in population trends before and after depletion (Table 4). Although the population experienced a considerable depression during the period of 2010-2012, it has since recovered (Figure 6a). We obtained similar results for ASIOTU, with an overall (1990-2021) stable population (S = 0.992 ± 0.010) and no significant difference between the periods (Table 4).
The average number of ASIOTU fledglings was 2.54 (95% bCI 2.38-2.72; n = 189) per successful nest. However, this population parameter declined over time (Figure 7c We did not find a correlation with the small mammal abundance in the year of breeding or the previous autumn ( Table 6).
The average number of BUBBUB fledglings was 2.26 (95% bCI 2.04-2.46; n = 81) per successful nest. This population parameter declined over time (Figure 7c Table 6). The effect of the Microtus vole abundance index in the previous autumn was also statistically significant and positive (Table 6). We had too few reliable observations of GLAPAS and AEGFUN breeding performance for analysis; therefore, we provide only a description of the average values: r = 2.75 (95% bCI 2.00-3.50; n = 8) and 1.75 (95% bCI 0.50-2.75; n = 4), respectively. Table 5. Description of owl breeding performance trends with small mammal cycles ("before") and since their depletion ("after").   On average, STRURA had 1.69 (95% bCI 1.58-1.80; n = 280) fledglings per successful nest. Despite the appearance of some differences in trends of breeding performance before and after small mammal cycle dampening (Figure 7b), they were not statistically significant ( Table 5). The overall trend of breeding performance was insignificant (β: −0.0014 ± 0.0122, t (24) = 0.112, p = 0.912; R 2 adj. = −0.0411; F (1;26) = 0.0125, p = 0.912). We did not find correlations with the small mammal abundance indices in the year of breeding or the previous autumn ( Table 6).

Small Mammal Monitoring
Overall, the small mammal densities and trapping indices in our study ( Figure 2) were similar to findings in the neighboring countries of Estonia [16,76] and Lithuania [77][78][79]. We found that meadow habitats had lower abundance of small mammals than forests (Tables A3 and A4). However, due to large within-class variation, no clear differences between age groups were found (Table A3). In Estonia [16] and Lithuania [78], increasing small mammal abundance has been recorded with increasing forest age in early meadowforest succession. Additionally, studies from Finland [80] and Norway [81] reported that mature forests had the highest abundance of voles.
However, a study conducted in northern Sweden suggested that young stands have higher small mammal diversity and abundance if a large amount of felling remains is left [82]. Many authors have found that high vegetation complexity, habitat structural diversity, and abundance of coarse woody debris are important factors that can help to ensure high diversity and abundance of small mammal species in young stands and unmanaged habitats under natural succession [79,[81][82][83][84][85].
The negative effects of intensive forestry have been found to be important at the landscape scale [83][84][85][86]. However, in a mosaic landscape, ecotones (with at least 100m buffer zone of habitat edges) have been found to contain the highest small mammal density and diversity [77,81,84].
Although the insufficient number of transects and occurrences in our study did not allow for statistical testing of forestry impacts, we found some declines in the small mammal numbers linked to forestry activities but unrelated to changes in other transects (Figure 2). Most of the small mammal monitoring transects were in intensively managed forests.
However, two of the longer-term areas were in protected areas (Apasalas and Žūklis) and also showed dampening of the cycles, suggesting larger than local (or management) effects on the dampening of the cycles.
One of the most robust explanations of cyclicity was provided by Hanksi et al. [17] with further extensions for different systems (see [87] for overview). One of those extensions, modelling multispecies rodent assemblages, revealed transient dynamics that alternated between long time periods with cyclic and non-cyclic fluctuations [19]. These fluctuations were expected to cover relatively small spatial scales, yet the phenomenon of dampened cycles was more recently found to occur Europe-wide [25], suggesting broader environmental drivers, for example, climate change [22,23]. However, in some parts of Europe, the period of dampened vole population cycles has been shorter than in others, refuting the generality of the climate forcing hypothesis [26]. Our results also showed clearly dampened vole cycles in Latvia (Figure 3a-c). The fact that such a dampening of cycles has not been reported in neighboring counties [14][15][16] suggests some smaller-scale processes, such as those described by predator-prey models. Although Hanski et al.'s models were created for the Fennoscandian environment, their generality has also been shown in central and western Europe [13]. According to these models and previous studies (see [87] for an overview), generalist predators tend to stabilize rodent dynamics, and nomadic avian predators have a similar effect on rodents, although they also increase the regional synchrony, whereas specialist predators have been thought to maintain the fairly regular rodent cycles [19].

Numerical Response of Owls
We found a statistically significant relationship between the proportion of voles (Microtus and bank voles) in owl diet and their relative abundance indices in nature (Table 3). This means that although we drew on mammal abundance information from autumn, it was still able to represent their abundance in owl prey. It is known that small mammal densities increase during summer [52] and that spring counts represent winter survival and reproduction [88], but the relative value of the current year (spring or autumn) still represents part of the cycle in cyclic environments [41]. We found preference for Microtus voles in every analyzed owl species in terms of the proportion in owl diet; these voles also accounted for a higher biomass proportion than bank voles (Table 3). Generally, a higher proportion of Microtus voles than bank voles in owl diet can be related to different breeding biology of voles and dispersal between vole species groups and predator-escaping behavior (see [19] for an overview).

Long-Eared Owl
We found ASIOTU to have the narrowest FNB among the investigated species. The calculated values were slightly lower than in Lithuania [28], possibly due to the pooling of the results to genus level. This species is known to be a small-mammal specialist in Europe [27] with a high proportion of Microtus voles in their diet [28,[89][90][91][92][93]. ASIOTU has been found to show strong functional responses of diet, breeding success, and dispersal to vole abundance [29][30][31]42]. It has even been suggested that species can adapt migratory behavior and breeding region selection during migration to account for vole abundance [94]. Moreover, this species may even exhibit repeated breeding attempts if the vole abundance is high [95,96].
Although the average breeding performance in Latvia was similar to the 2.94 ± 0.42 (µ ± SD; n = 1339) recorded in Finland [97] and 2.39 (n = 72) in the United Kingdom [95], we observed a significant decline in the number of fledglings per successful nest, i.e., more than one chick in three generations (5.7 years; [98]). The steepest decline occurred in the last two generations and matched the time of dampened populations of small mammals (Figure 7c). The declining breeding performance did not have an impact on the population change (2004-2016), but extension of the trend (2004-2021) showed a significant decline (β: −0.0530 ± 0.0191, t (16) = −2.772, p = 0.0136). We suspect that for a longer period of time, the ASIOTU population was supported by immigration of migrating individuals hatched elsewhere [29,90,99,100] and that a later decline implies a delayed response of returning individuals of Latvian origin.
Habitat composition and prey abundance have been found to be the most important factors shaping local ASIOTU populations [90,[101][102][103] because the species shows no strong territorial defence and hunting grounds may largely overlap between neighboring pairs [104,105]. We found no correlation between the breeding performance and prey abundance indices in the year of breeding or the previous year (Table 6), likely due to selection of breeding territories with sufficient abundance of prey. This is supported by the knowledge of species benefitting from relatively small landscape elements, for example, flower strips [103].

Tengmalm's Owl
Although we had only two samples of AEGFUN diet, its FNB suggested a high specialization, which was supported by a high proportion of voles in the diet. The observed proportions in Europe have shown a high importance of voles (overview in [32]), averaging 54.89% according to previous studies. This species has a strong functional response to vole abundance influencing habitat selection via hunting behavior [47,[106][107][108], the timing of breeding and breeding performance [6,39], and survival [39,43]. Even with a certain degree of carry-over effect, the species has shown strong adaptability to fluctuating food conditions in terms of breeding performance [38].
The average long-term number of fledglings in Finland is 4.04 ± 0.62 (µ ± SD; n = 13,817) [38] and around 2 fledglings per successful nest in poor vole years [32,109]. The scarcely available data on the breeding performance in Latvia suggests that it is similar to that in vole depression years in Finland.
We found a steep decline of AEGFUN population throughout the studied period (Figure 6d), but it was slower with depleted population dynamics of small mammals (Table 4). We expected this to be an artefact of some better seasons or immigration from Fennoscandia and Russia [32] rather than an actual difference; therefore, we extended the period of analysis in the "after" group. Our results (β: −0.0618 ± 0.0074, t (16) = −8.302, p = <0.0001) showed a decline since 2004 closely matching the overall population decline and the slope of the period with pronounced population dynamics of small mammals. Some researchers have hypothesized potentially negative effects of increasing STRURA population on the AEGFUN population [32]. However, we did not find any AEGFUN as a prey of STRURA, although superpredation is known [34] and both species coexisted in the same study areas (authors' personal observations). Furthermore, in central Europe, breeding in the proximity of STRURA has been found to protect AEGFUN against STRALU [110][111][112].
Population declines have also been reported in Finland, Sweden, and Estonia [113], suggesting larger-scale factors affecting the population. This species is a mature spruce and mixed forest specialist [47,101,108,[114][115][116][117][118][119]. These are habitats with some of the highest densities of small mammals [79][80][81][82][83][84][85]. We consider the loss of species-specific habitats to be the most important factor in population decline, amplified by dampened dynamics of small mammal populations in Latvia. The forestry intensity, as measured by tree cover loss, is increasing in Latvia and, in particular, in priority sites for species conservation [120].

Eurasian Pygmy Owl
We found average level of specialization of GLAPAS and the preference for Microtus voles was strongest among the analyzed owl species (Table 3), although with a low proportion of voles in the diet. The vole proportion was similar to the breeding season diet in Finland [34] and in central Europe [121]. Masoero et al. found a strong numeric and functional response of GLAPAS to vole abundance in winter [33], suggesting not only ageand gender-specific preference for voles but also stronger migratory behavior during lowvole-density years in boreal Finland. During the years of higher vole population densities, breeding density and performance of GLAPAS also increase [6,122]. The dependency on voles has been found to be stronger in boreal than boreonemoral regions, with breeding both in low and peak vole years in the latter [40]. In the boreonemoral zone, the onset of breeding was later with no correlation with breeding performance, and the clutches were slightly smaller than in the boreal zone [40].
For the few documented records of successful breeding in Latvia, the values were markedly lower than 5.85 ± 0.55 (µ ± SD; n = 13,817) in Finland [97] and boreal Norway (6.9 ± 1.1) and somewhat lower than in boreonemoral Norway in vole crash years (3.7 ± 2.8) [40]. The difference relative to boreonemoral Norway indicates a possible cumulative effect of longer-term dampened population cycles of small mammals, which is supported by a declining GLAPAS population. The population in Estonia and Lithuania is increasing [113], but it is declining in Latvia (Figure 6e and Table 4) and in Finland [97,113]. It has to be noted that only Finland and Latvia were able to provide analytical assessment of the population (type: interval) in the last report on the Article 12 of the Birds Directive [113]; therefore, it cannot be ruled out, that the increase in the other Baltic states is more based on increased survey efforts and knowledge than a genuine change. Although irruptions linked with low rodent availability occur from time to time [123], it is unlikely that Finland and Latvia are source populations for neighboring countries with declining populations themselves, despite the increased distribution of the species [51]. This is supported by similar patterns of yearly indices in Latvia ( Figure 6e) and Finland [97] but with a steeper decline in Latvia.
GLAPAS is known to be a structurally rich, mature spruce and mixed forest specialist species during the breeding season [122,[124][125][126][127][128][129][130][131], and clearcuts and logging have been shown to affect habitat suitability [131], as well as population size [132]. Structurally rich mature forests are habitats with some of the highest densities of small mammals [79][80][81][82][83][84][85]133]. Latvia and Finland are the countries in Europe with the highest forestry activity, even in protected areas [134]. We suspect the loss of species-specific habitats to be the most important factor in the population decline, amplified by the dampened dynamics of small mammal populations in Latvia. The forestry intensity, as measured by tree cover loss, is increasing in Latvia and, particularly in priority sites for species conservation [120].

Ural Owl
One of the highest and temporarily increasing FNB values was found in STRURA. We found a low proportion of voles in the diet of this species. This proportion, when compared by count, was lower than in Finland [34,35,135,136], similar to that in Belarus [137][138][139], and higher than that in Slovenia [35]. When comparing the food niche as a whole, STRURA diet in Latvia was found to be similar to that in Finland during the low vole phase [35]. Although this species is known to be a generalist predator, a strong functional response to vole abundance has been proven in Fennoscandia, ranging from the timing of breeding and breeding performance [3,6,39,41,140] to winter survival [3,6,39,41,45] and even demonstrating a carry-over effect from the previous year [48] and a change in behavior [46,141,142].
Not only the food niche but also the breeding performance of STRURA in Latvia was similar to that in Finland in bad vole years. In Latvia, we observed, on average, 1.69 (95% bCI, 1.58-1.80; n = 280) fledglings per successful nest and no temporal trend. The corresponding overall value in Finland (1986-2016) is 2.59 (±0.43 SD, n = 18,901; [98]) and between 1.3 and 2 [3,41] in bad vole years, roughly matching our results. Given the strong numerical response to voles, we expected a declining trend in breeding performance, but we did not find it. We consider this an example of strong parental investment [142][143][144] as reflected by adjustments in hunting activity and possibly habitat selection [101], demonstrating the high plasticity of the species [35]. As Figure 2 shows, even with dampened small mammal cycles, there are habitats with high prey abundance, allowing prey to meet the demands of the young. The size of nest boxes in Latvia is similar to that in Finland [35] and cannot be suspected as a reason for lower breeding performance.
Increasing STRURA populations and expanding range, even increasing the niche of utilized habitats, was observed in many parts of Europe during the first decade of the 21st century [145][146][147][148]. This overlaps with the increase in Latvia and breeding occurrences in a mosaic landscape [101]. Given the extent of population increase, some unknown large-scale factors are most likely the explanation. Nevertheless, in Latvia, the period of steep decline in the species population overlapped with the dampening of the small mammal population dynamics. We consider the relative abundance of small mammals to be an important collider for overall habitat change, as species ecological niche analysis in Latvia suggests strong dependency on large forest massifs with dominance of mature forests and only some openings [102]. These are habitats found to be important for the same species elsewhere [146,149,150]. Although the range is still expanding in Latvia [151], the overall population size is declining [113]. The forestry intensity, as measured by tree cover loss, is increasing in Latvia, in particular in priority sites for species conservation [120]. We consider this as an argument for the conservation of mature forests important for the species and holding higher densities of its main prey, i.e., small mammals [81].

Tawny Owl
The highest FNB value was found in STRALU with a relatively low proportion of voles in the diet. The average FNB value was slightly lower than in Lithuania [28]. We observed a temporal increase in FNB, which was similar to the observation in Lithuania, with a declining proportion of Microtus voles [36]. This species is known to be a generalist [152]. Its food composition can considerably vary between breeding regions within the same year and between years in the same breeding territory [28,36,[152][153][154][155][156][157][158][159][160]. Nevertheless, in the cyclic environment of Fennoscandia, a strong numerical response to vole abundance has been reported, including the timing of breeding [6,160], breeding performance [3,5,6,39,160,161], and winter survival [3,39].
Both the population change and breeding performance of STRALU were stable and showed no differences relative to pronounced and dampened vole cycles. Breeding performance was lower than 3.26 ± 0.41 (µ ± SD; n = 9668) in Finland, where the population was also stable [97], as well as in Lithuania, where an increasing trend of breeding performance (2002-2014) was observed and co-occurred with a decline in the number of breeding pairs [36]. We consider the relatively low breeding performance in Latvia to be related to the high population density (estimated to be around 16,604 in Latvia and below 4000 in Lithuania) [113]. The observed depression of the STRALU population from 2010 to 2012 in Latvia partially matched with Lithuania [36]. This was likely the consequence of two consecutive snow-rich winters, with multiple freeze-thaw events forming ice sheets in snow cover-factors reducing species survival [45,162,163]. This event did not affect breeding performance, and the population recovered quickly.
We think that the quick population recovery and overall stable breeding performance, even with increasing FNB values, was possibly due to the breeding habitat availability. Although this species is a well known generalist breeding from cities to large forest massifs in more southern latitudes [34], in the boreonemoral region, STRALU has been found to prefer forest edges over the interior [101,164]. With increasing forestry and forest fragmentation, more suitable landscapes for this species are created [101], probably overwhelming the negative effects of depleted small mammal populations.

Eagle Owl
The largest European owl species, BUBBUB, is known to be a generalist predator, with the proportion of rodents in its diet ranging from 0 to 97.7%, with average of 49.7% among 182 studies (overview in [37]). We do not have reliable information on the diet of BUBBUB in Latvia, but during the ringing of the young, many bird feathers are found, as most of the known breeding sites are in close proximity to waterfowl lakes and landfills [101]. This species is resident with no known seasonal migrations in Europe [37,165], and breeding dispersal occurs mostly due to the loss of a mate [37]. We speculate that the BUBBUB Life 2023, 13, 572 20 of 28 population in Latvia depends highly on voles, at least in winter, when bodies of open water are typically frozen and most waterfowl and gulls in Latvia have moved to wintering areas [151]. This is supported by a study conducted in Finland evaluating the robustness of the alternative prey hypothesis for BUBBUB [166]. A correlation was found between vole abundance in nature and their proportion in the diet, and the proportion of alternative prey was found to be nearly independent of its abundance in the field [166]. Several other studies highlight the high proportion of voles in the BUBBUB diet [167][168][169].
The overall average number of fledglings per successful nest in Latvia was similar to that in Europe-around two (overview in [37]); however, we observed a declining trend, with a loss of more than one fledgling in two generations (generation length is 12.1 years [170]). We found that breeding performance was correlated with the abundance index of small mammals in field, and the correlation with the value from the previous autumn was stronger (Table 6). Other studies found that BUBBUB pairs with a diet based on high-value foods (rabbits and rodents) have comparatively larger broods and breed earlier [171,172], and higher reproductive productivity was associated with a higher proportion of the main prey (rats and rabbits in Spain) in the diet [173]. Our data do not lend themselves to examining such a relationship between breeding performance and diet. Nevertheless, we assume that the negative effect of a reduced abundance of small mammals highlights a carry-over effect, influencing adult fitness in spring and thereby reducing the breeding performance. This phenomenon is well known in STRURA [41,48,143,174,175] and has proven to be of increasing importance with the size of an owl species in Finland [6]. Ecological niche analysis in Latvia also suggests the importance of habitats with higher vole abundance [81] for BUBBUB [101].
Although we had annual data on a limited number of nests, they formed an important part of the whole population estimated at around 24 breeding territories, indicating a declining national population trend [113]. We consider our findings of a possible carry-over effect to be important in species conservation and to be linked to population decline via reduced breeding performance, as well as reduced winter survival, as both should be related via fitness, although this relationship needs to be studied more directly. Nevertheless, we consider the conservation of habitats important for breeding and winter feeding, together with nest site protection from ground predators, which is necessary to reduce the effects of dampened population dynamics of small mammals.

1.
Small mammal relative abundance indices have shown depleted population cycles since approx. 2004. This has impacted the breeding performance, food niche breadth, and population trends of owl species to various degrees depending on the particular species; 2.
The number of ASIOTU fledglings has declined since the depletion of small mammal populations. The population size of the species declined later and was significant for the period from 2004 to 2021. ASIOTU is the most specialized of the analyzed owl species in terms of the proportion of voles in the diet; 3.
The breeding performance of the three forest specialist species AEGFUN, GLAPAS, and STRURA in Latvia was similar to vole depression years in the boreal and boreonemoral regions; 4.
Populations of GLAPAS and AEGFUN declined in Latvia and showed no difference compared to periods with pronounced or depleted population dynamics of small mammals. In contrast, the population of STRURA has shown a significant decline since rodent depression. We suggest the depletion of the small mammal population dynamics to be an important negative contributing factor to more important effects of forestry, although the impact of forestry needs to be investigated further; 5.
Neither the breeding performance nor population size of STRALU changed between the compared periods with pronounced and depleted population dynamics of small mammals. This suggests a strong plasticity of the species, as food niche breadth was temporarily increased; 6.
We found evidence that suggests the dependency of BUBBUB on voles via a carry-over effect. The breeding performance of BUBBUB was significantly correlated with the abundance indices of small mammals in nature in the previous autumn.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/life13020572/s1, Table S1: Description of owl diet per year; Table S2: Population indices of small mammals; Table S3: Population indices of owls; Table S4: Description of owl breeding performance.